Route of hydrodynamic turbulence in rotating shear flow: Secondary perturbation 
triggers growing eigenmodes through the elliptical instability 
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Origin of hydrodynamic turbulence in rotating shear flow, e.g. plane Couette flow including the 
Coriolis force, is a big puzzle. While the flow often exhibits turbulence in laboratory experiments, 
according to the linear perturbation theory it should always be stable. We demonstrate that the 
secondary disturbance to the primarily perturbed flow triggers the elliptical instability in such a 
system and hence an exponential growing eigenmode. This result has an immediate application 
to astrophysics and geophysics which often exhibit a turbulent flow with the Keplerian angular 
momentum profile. We address the origin of turbulence in such a Keplerian flow which is similar to 
rotating Couette flow. 
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Despite of many efforts devoted, the origin of hydro- 
dynamic turbulence is still poorly understood. One of 
the main problems behind this is that there is a sig- 
nificant mismatch between the predictions of linear the- 
ory and experimental data. For example, in the case of 
plane Couette flow, laboratory experiments and numeri- 
cal simulations show that the flow may be turbulent for 
a Reynolds number, R ~ 350, while according to the lin- 
ear theory the flow should be stable for all R. Similar 
mismatch between theoretical results and observations is 
found in geophysical and astrophysical contexts, where 
the accretion flow of neutral gas with Keplerian angular 
momentum profile is a common subject. A Keplerian ac- 
cretion disk flow having a very low molecular viscosity 
must generate turbulence and successively diffusive vis- 
cosity, which support transfer of mass inward and angu- 
lar momentum outward. However, theoretically this flow 
never exhibits any unstable mode which could trigger 
turbulence in the system. On the other hand, laboratory 
experiments of Taylor-Couette systems, which are simi- 
lar to Keplerian disks, seem to indicate that although the 
Coriolis force delays the onset of turbulence, the flow is 
ultimately unstable to turbulence for Reynolds numbers 
larger than a few thousand Q], even for subcritical sys- 
tems. Longaretti reviews the experimental evidence 
for the existence of turbulence in subcritical laboratory 
systems, and concludes, based on phenomenological anal- 
ogy, that a similar process must happen in astrophysical 
accretion flows. Indeed, Bech & Anderson Q see tur- 
bulence persisting in numerical simulations of subcritical 
rotating flows for large enough Reynolds numbers. 

How does shearing flow that is linearly stable to per- 
turbations switch to a turbulent state? Since last decade, 
many authors have come forward with a possible ex- 
planation of this fact, based on bypass transition (see 
IE IE IE IE and references therein) , where the de- 
caying linear modes show an arbitrarily large transient 
energy growth at a suitably tuned perturbation. In lieu 
of linear instabilities e.g. magnetorotational instability, 



the transient energy growth, supplemented by a non- 
linear feedback process to repopulate the growing dis- 
turbance, could plausibly sustain turbulence for large 
enough Reynolds numbers. 

However, in the case of rotating shear flow, e.g. 
Couette- Taylor flow, Keplerian accretion flow, the tran- 
sient energy growth is insignificant for three-dimensional 
perturbations. The Coriolis effect, which absorbs the 
pressure fluctuation, is the main culprit to kill any growth 
of energy. Nevertheless, at a very large Reynolds num- 
ber, certain two-dimensional perturbations may produce 
large transient growth in such a flow. However, in 
two-dimension, the underlying perturbations must ul- 
timately decline to zero in the presence of viscosity 
[E 0| • To overcome this limitation, it is necessary to 
invoke three-dimensional effects. Various kinds of sec- 
ondary instabilities, such as the elliptical instability, are 
widely discussed as a possible route to self-sustained 
turbulence in linearly perturbed shear flows (see, e.g. 
EEmHHQIlSlIEIH E3) These effects, 
which generate three-dimensional instabilities of a two- 
dimensional flow with elliptical streamlines, have been 
proposed as a generic mechanism for the breakdown 
of many two-dimensional high Reynolds number flows 
whose vortex structures can be locally seen as ellipti- 
cal streamlines. However, to our knowledge, such effects 
have not been discussed properly in literatures for ro- 
tating Keplerian flows which have a vast application to 
various natural phenomena. More than a decade ago, 
the possible role of the elliptical instability in astrophys- 
ical accretion disk physics was first explored [19( and 
it was shown that angular momentum may be trans- 
ferred from the disk to the tidal source by the instability 
effect [IE EI- Therefore, naturally we are motivated 
to see whether these three-dimensional instabilities are 
present in rotating shear flows which consist of ellipti- 
cal streamlines under two-dimensional perturbation. We 
plan to show that in presence of secondary effects, three- 
dimensional perturbation can generate large growth and 
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presumably trigger turbulence in rotating shear flows. 
Unlike the growth of transient kind which is significant 
in two-dimension only, in the present case it is essen- 
tially the three-dimensional effects which produce large 
growth. Possibility of large growth in shear flows with 
rotation by three-dimensional perturbation opens a new 
window to explain hydrodynamic turbulence. 

Let us consider a small patch of rotating shear flow 
whose unperturbed velocity profile corresponds to a lin- 
ear shear of the form Uo = (0, — x, 0). Because of rota- 
tion, a Coriolis acceleration acts on the fluid element hav- 
ing angular frequency il — (0,0, f^) where according to 
choice of our units £^3 = 1/q. The parameter q is positive, 
describing the radial dependence of Q(r) = Qo(r Q /r) q 
(see for detail). We choose all the variables to be di- 
mensionless. When the Reynolds number is very large, 
this flow exhibits a large transient growth under a suit- 
able two-dimensional perturbation [9( which modifies the 
linear shear profile as 



U = Uq + u? = (w x , — x + w y , Q) = A.d, 



(1) 



fc k 
w x = C-f sm(k x x + k y y), w y — — C~f sin(k x x + k y y), 



where d is the position vector and A is a tensor of rank 
2. Now expanding the sin(fc x x + k y y) terms and choosing 
a small patch, A is given by 



A = 4= I -(1 + Ce) 



C( l - g) 




(2) 



where e = {k x / k) , k = yk 2 + k y , and C is the ampli- 
tude of vorticity perturbation. Here k x = k X Q+k y t, which 
basically is the radial component of wave-vector, varies 
from —00 to +00, where k x o is a large negative number. 
The above plane wave typed perturbation is frozen into 
the fluid and is sheared along with the background flow. 
At t = 0, the effective wave vector of perturbation in the 
x direction (k x ) is negative, which provides very asym- 
metric leading waves. As time goes on, the wavefronts 
are straightened out by the shear and \k x \ decreases. At 
the time when the transient growth is maximum, k x ~ 
and the wavefronts become almost radial. At yet later 
time, the growth decreases and the wave becomes of a 
trailing pattern. Clearly U describes a flow having gen- 
eralized elliptical streamlines with e, a parameter related 
to the measure of eccentricity (2^|, runs from to 1 as 
the perturbation evolves. Now we plan to study how 
does this perturbed shear flow behave under the further 
perturbation, namely a secondary perturbation. 



The linearized equations for the evolution of a sec- 
ondary perturbation w, such that U — > U + u, to the 
flow are 



(d t + U.V^j u + u.VU + 2nxu = - Vp + V z u, 

V.u = (3) 
where p is pressure in the fluid element and R is the 
Reynolds number. A secondary perturbation of the plane 
wave kind is given by 
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(Ui,p) = (vi(t),p(t))exp(ik m (t)x m ), 



(4) 



where the Latin indices, run from 1 to 3, indicate the 
spatial components of a variable, e.g. x m = (x,y,z). 
Therefore, after some algebra, we obtain the evolution of 
a linearized secondary perturbation 



ij + A) v k + 2 e mk3 n m v k = ~ip kj - fc 2 , (5) 



k n V r , 



(6) 



where 'over-dot' indicates a derivative with respect to t, 
tmkj is a Levi-Civita tensor, and k 2 = k m k m . A simi- 
lar set of equations was obtained by Bayly except 
that they now have additional terms induced due to the 
Coriolis and the viscous effects along with a modified A. 
For a fixed k x , the components k\ and k^ of wave- vector 
[fc = k m = (ki, &2) ^3)] of secondary perturbation oscil- 
late in time with the angular frequency w = ^C(l — e )i 
while the vertical component, fc3, is constant. As we 
choose the signature of the background Minkowski space- 
time to be [—,+,+,+], it does not matter whether any 
Latin index appears as a lower case or an upper case. For 
example, Aj — Ajk, where j and fc indicate the row num- 
ber and the column number of the associated matrix re- 
spectively. Projecting out eqn. JSJ) by P? = 5? — k~ 2 W ki 
and using eqn. JSJ), we obtain 



[2e mkj il v + — fc 



R ) fc 2 



(7) 



Now we specifically concentrate on the flow having low 
viscosity. Therefore, we neglect the viscous term in eqn. 
(Q comparing others and rewrite 
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A,-i = 



2k 1 k 2 



An + 



- lj A n + 

2k~ 



2k 1 k 2 
~k^~ 



^A 11 + 



(A 2 i + Sl 3 ) 
)- 

— (A 2 , + Q ;i ) 



(A 2 i + Sl 3 ) - A21 - 2f2 3 

2fc 2 fe3_ 



r 



^ (a 12 - n 3 ) + 2^3 - A12 + 

(A 12 -n 3 ) + (^-i) a 22 

, 2koki a 

22 



~k^~ 



(Avt - fi 8 ) + 



22 \ 



0/ 



(8) 



The general solution of eqn. |JSJ for a particular can 
be written as a linear superposition of Floquet modes 



Vi(t) = exp(at) fi (</>), 



(9) 



where <p — vot, fi(4>) is a periodic function having time- 
period T = 2Tt/m 1 and ct is the Floquet exponent which 
is the eigenvalue of the problem. As k x varies, e changes, 
therefore the eigenvalue changes. Clearly, if a is positive 
then the system is unstable. To compute a, one has to 
evaluate the elements of a matrix M m i(2ir) whose eigen- 
values and eigenvectors are respectively e aT and /,-(27r), 
while M m i((j>) itself satisfies 



(10) 



with the condition Mjj(0) = Sji. Clearly, l|10|l is the evo- 
lution equation of velocity. If /1 is a real positive eigen- 
value of M m i(2it), then the corresponding a — ln(fj,)/T. 
One can evaluate fi by an elementary numerical tech- 
nique. However, eqn. i jlOjl has exact analytical solution 
for an initial perturbation ko = (0,0, 1). In this case, k 
remains constant throughout for a fixed e and therefore 
the growth rate, a, is the highest eigenvalue of A™ given 
by 



<r = ^£"(2^3-1) (2ft 3 -0- 



(11) 



• When fl 3 = 0, a = y/C(e - 1). This verifies that 
non-rotating two-dimensional shear flow is always hydro- 
dynamically stable under a pure vertical perturbation. 

• When SI3 = 1/2, cr = \/(e. Therefore, rotating shear 
flow with a constant angular momentum profile is always 
hydrodynamically unstable. The energy growth rate of 
perturbation increases with the strain rate, i.e. eccen- 
tricity, of the flow. 

• When fl 3 = 2/3, a = y/(e - (4 - 3<)/9. Therefore, 
a Kcplcrian flow with elliptical streamlines is hydrody- 
namically unstable under a pure vertical perturbation, 
only if C > 1/3. 

However, apart from the vertical one, there are some 
other three-dimensional perturbations |23| which can 
generate instability in rotating shear flows with £ < 1/3, 
that we describe by numerical solutions. We essentially 
focus on the Keplerian flow when q = 3/2. As primary 
perturbation evolves with time, eccentricity decreases, 
and then energy growth rate due to secondary pertur- 
bation changes. Figure^ shows the variation of max- 
imum growth rate, cr max , as a function of eccentricity 
parameter, e, for various values of £, the amplitude of 
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FIG. 1: (a) Variation of maximum growth rate as a function 
of eccentricity parameter when solid, dotted, dashed and long- 
dashed curves indicate the results for £ = 0.2,0.1,0.05,0.01 
respectively, (b) Variation of optimum growth rate as a func- 
tion of vertical component of perturbation wave-vector when 
various curves are same as of (a), q = 3/2 throughout. 



primary perturbation, in the Keplerian flows. By "maxi- 
mum" we refer the quantity obtained by maximizing over 
the vertical component of wave- vector, £3. While a ver- 
tical perturbation gives rise to the best growth rate for 
C > 1/3, as we show above analytically, for £ < 1/3 
other three-dimensional perturbations maximize growth 
rate. At small e and large (, the streamlines of the flow 
essentially become circular (see eqn. (j2J), and thus in 
absence of any significant elliptical effect growth rate 
severely decreases. On the other hand, when e and £ 
both are small, the background reduces to plane shear 
structure and therefore any growth arises due to primary 
perturbation only. 

Figure^) shows the variation of optimum growth rate, 
cr op t, as a function of k 3 for various values of £, in the 
Keplerian flows. By "optimum" we refer the quantity 
obtained by maximizing over e. Interesting fact to note 
is that the optimum growth rate is obtained always for 
three-dimensional perturbation with significant vertical 
effect, i.e. a non-zero value of k 3 . Moreover, as £ in- 
creases, instability occurs at larger e with higher £3. 
Therefore, three-dimensional elliptical instability is more 
prompt at larger £. 

Above results verify that at a range of e the instan- 
taneous three-dimensional growth rate due to secondary 
perturbation is always real and positive, which motivates 
us to compute the actual growth itself during simulta- 
neous evolution of both the perturbations. As primary 
perturbation evolves, k x varies with time and therefore 
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A does so. Thus, in reality, the wave- vector of secondary 
disturbance is no longer periodic, though the correspond- 
ing equation JTjJ) is still valid. To compute growth in en- 
ergy one has to find out the elements of a matrix 



M tk (t) = M im (t) T M mk (t) 



(12) 



whose largest eigenvalue is growth in energy at a time t. 
Clearly, M.%k(f) is the instantaneous energy of the flow. 
k x and k\ start with a large negative value when the 
flow is highly eccentric. With time k x (as well as ki) 
decreases in magnitude and finally becomes zero when 
the streamlines become circular. If £ = 0, then k x is 
same as k%. However, for £ > 0, k\ increases faster than 
k x , as follows from eqn. ijBJ), and becomes positive during 
the evolution of perturbation. When k\ — > 0, growth 
maximizes. 

Figure shows the variation of growth as a func- 
tion of t for various values of (. Note that as £ in- 
creases, growth maximizes at an earlier time. Figure 
shows the variation of maximum growth as a function 
k-sjki- The quantity kj,jki carries information about 
how three-dimensional the flow is. We know that in 
two-dimensional plane shear flow the maximum growth, 
G m ax, scales with fc 2 0. However, for ( > 0, G max 
decreases at small k^/k2, while increases at large fea/fo- 
This clearly proves that three-dimensional secondary per- 
turbation triggers elliptical instability which produces 
larger growth. As k x0 and k w (~ i? 1 / 3 0,0) increase, 
the effects due to elliptical instability increase, and thus 
the corresponding growth does so. When k x o — kio — 
-10 3 , G max - 4 x 10 4 at k 3 /k 2 = 1 for ( = 0.1, which 
is an order of magnitude larger compared to that for 
C = 0. If we consider a smaller R with the corresponding 
k x o = fcio = — 10 2 , then G max at k^lk^ = 1 decreases 
to ~ 2 x 10 3 for £ = 0.1, which is still a factor of two 
larger compared to that for £ = 0. Therefore, this is 
confirmed without any doubt that three-dimensional el- 
liptical instability efficiently triggers turbulence in shear 
flows. 

When t increases from to t max = —k x o/k y , k x (t) 
varies from k X Q (~ — oo) to 0, and e decreases from 
1 to 0. Note that G max ~ 2 x 10 3 - 4 x 10 4 when 
kf ~ R ~ 10 6 — 10 9 . As this large growth is the re- 
sult of three-dimensional perturbation, even in presence 
of viscosity the underlying perturbation effect should sur- 
vive. Presumably, this growth factor is enough to trigger 
non-linear effects and turbulence in such flows. There are 
many important natural phenomena where the Reynolds 
number is very large. In astrophysical accretion disks R 
always could be 10 10 — 10 14 |9( because of their very low 
molecular viscosity. Therefore, the present mechanism 
can have a very good application to such disk flows to 
explain their turbulence puzzle when it is cold and neu- 
tral in charge. On the other hand, we argue that the 
subcritical transition to turbulence in Couette flow may 
be the result of secondary perturbation which triggers 
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FIG. 2: (a) Variation of growth as a function of time due 
to simultaneous evolution of both perturbations, when &3 = 
0.8. Solid, dotted, dashed and long-dashed curves indicate the 
results for £ = 0,0.01,0.05,0.1 respectively, (b) Variation of 
maximum growth as a function of k^/k2, when various curves 
are same as of (a). Other parameters are k x o — &io = —1000, 
feo = 1, and q = 3/2. 



elliptical instability modes in the system. It is to be 
seen now whether all shear flows, exhibit subcritical tur- 
bulence in the laboratory, do also exhibit large growth 
due to secondary perturbation, according to the present 
mechanism. 

The author is grateful to Ramesh Narayan for suggest- 
ing this problem and for extensive discussion and encour- 
agement throughout the course of the work. This work 
was supported in part by NASA grant NNG04GL38G 
and NSF grant AST 0307433. 



[10 

[11 
[12 

[is; 

[14 
[15 
[16 



D. Richard, & J.-R Zahn, A&A 347, 734 (1999). 
R Longaretti, ApJ 576, 587 (2002). 

K. Bech, & H. Andersson, J. Fluid Mech. 347, 289 
(1997). 

K. Butler, & B. Farrell, Phys. Fluids A 4(8), 1637 (1992). 
S. Reddy, & D. Henningson, J. Fluids Mech. 252, 209 
(1993). 

L. Trefethen, A. Trefethen, S. Reddy, & T. Driscoll, Sci- 
ence 261, 578 (1993). 

G. Chagelishvili, J. -P. Zahn, A. Tevzadze, & J. Lomi- 

nadze, A&A 402, 401 (2003). 

O. Umurhan, & O. Regev, A&A 427, 855 (2004). 

B. Mukhopadhyay, N. Afshordi, & R. Narayan, ApJ 

629, 383 (2005); N. Afshordi, B. Mukhopadhyay, & R. 

Narayan, ApJ 629, 373 (2005). 

R. Pierrehumbert, Phys. Rev. Lett. 57, 2157 (1986). 

B. Bayly, Phys. Rev. Lett. 57, 2160 (1986). 

A. Craik, & W. Criminale, Proc. R. Soc. London Ser. A 
406, 13 (1986). 

M. Landman, & P. Saffman, Phys. Fluids 30(8), 2339 
(1987). 

C. Hellberg, & S. Orszag, Phys. Fluids 31(1), 6 (1988). 
F. Waleffe, Phys. Fluids A 2(1), 76 (1989). 

A. Craik, J. Fluid Mech. 198, 275 (1989). 



■5 



[17] S. Le Dizes, M. Rossi, & K. Moffatt, Phys. Fluids 8(8), 
2084 (1996). 

[18] R. Kerswell, Ann. Rev. Fluid Mech. 34, 83 (2002). 
[19] J. Goodman, ApJ 406, 596 (1993). 

[20] S. Lubow, J. Pringle, & R. Kerswell, ApJ 419, 758 
(1993). 

[21] D. Ryu, & J. Goodman, ApJ 422, 269 (1994). 

[22] Note that e is a parameter related to the measure of ec- 



centricity but not the eccentricity itself. 
[23] By vertical perturbation we mean that only the ver- 
tical component of initial perturbation wave-vector is 
non-zero, while any perturbation with a non-zero ver- 
tical component of initial wave-vector is called three- 
dimensional perturbation. 



